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We review the consequences of the growth and evolution of Black Holes on the distribution of stars 
and Dark Matter (DM) around them. We focus in particular on Supermassive and Intermediate- 
Mass Black Holes, and discuss under what circumstances they can lead to significant overdensities 
in the surrounding distribution of DM, thus effectively acting as DM Annihilation Boosters. 
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•rH . I. INTRODUCTION 

^ • 

■ To identify the nature of Dark Matter (DM) particles, three different experimental strategies have been devised. 
First, high-energy colliders like the Large Hadron Collider at CERN pj may soon be able to produce DM particles, 
and possibly identify them (see e.g. Refs. 0] and Q and references therein). Another possibility is Direct Detection, 
based on the search for rare events in which a DM particle scatters off the nuclei of large detectors (see e.g. Ref. 

for a review). Finally, Indirect Detection is based on the search for the products of annihilation or decay of DM 
particles (such as photons, neutrinos and anti-protons). Many candidates have been proposed for DM, but the most 
studied ones belong to the class of Weakly Interacting Massive Particles (WIMPs) [1, Q , that includes, among others, 
the lightest neutralino in Minimal or Non-Minimal Supersymmetric extensions of the Standard Model , and the 
lightest Kaluka-Klein state in models with Unified Extra Dimensions [1, [l3| • 

Here we focus on Indirect DM searches. Ideal targets for Indirect DM searches are nearby, high DM density 
regions, such as local celestial bodies like the Earth and the Sun, where DM can be captured, or centers of galactic 
halos (see e.g. Refs. [ill El! 113, [IE [Hj [I3|)- Two classes of astrophysical objects of particular interest are 
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SuperMassive Black Holes (SMBHs) [Tl,[Ti, with masses from 10^ to IO^Mq and the more speculative Intermediate 
Mass Black Holes (IMBHs), with a mass from 20Mq to 10^ Mq (see e.g. Refs. [13, [H], [Hi and references therein). 
Both these classes of compact objects can influence the distribution of DM in which they are embedded, leading to 
strong overdensities: we review here the impact of the formation and growth of BHs on the surrounding distribution 
of matter, and the consequences for Indirect DM searches. Being the annihilation flux proportional to the integral of 
the DM density squared, scenarios where the density is boosted by the presence, or the growth, of a central BH, are 
very promising for Indirect searches. BHs can thus be considered as DM Annihilation Boosters. 

The paper is organized as follows: in Section|TT]we discuss our current understanding of DM proflles in halos without 
central BHs. In Section lllll we address the case of a DM distributon evolving around an already-formed BH. In Section 
IIV[ we discuss the case of adiabatic growth and the conditions under which DM "spikes" can form. In Section IVl we 
will focus on IMBHs, and show that some of the problems associated to the formation of spikes around SMBHs can 
be evaded, with important consequences for Indirect DM searches. Conclusions will be drawn in Section VI. 



II. DARK MATTER PROFILES WITHOUT BLACK HOLES 



There is strong evidence in favor of the presence of SMBHs at the center of every galaxy with a substantial bulge 
component [H, |T^, '23^1 and it has been suggested that even globular clusters can harbor IMBHs . Throughout 
the paper, we will generically refer to large gravitationally bound systems like globular clusters, galaxies and clusters 
of galaxies, as galaxies, and we will refer to their central region as the nucleus, which may host BHs. Such compact 
objects account roughly for 10~'^% of the baryonic mass of the galaxy, which is composed of stars, intergalactic dust 
and DM. We focus on a generic WIMP scenario, with an annihilation cross section of order av w 10~^^cm^s~^ and a 
mass ranging between the GeV and the TeV scale. 

Since we want to characterize how BHs influence the surrounding distribution of matter, we need to specify how 
DM is distributed before the BHs form, and use this information as initial condition for the problem at hand. Profiles 
without any central object also receive particular attention per se, since the cuspiness of a DM halo without BH can 
give informations about the "coldness" of the DM candidate 24]. Recently, A^-body simulations of galaxies found 
results in favor of power law profiles (with slope from —1 to —1.5) for the nuclear region, emphasizing the contrast 
with direct observations, such as rotation curves of Low Surface Brightness galaxies (LSBs) [2^ [2^ [23| and X-ray 
imaging, which s ugg est instead the presence of flat DM cores. 

Navarro et al. [28| and Reed et al. [2§| used A''-body techniques to simulate the high resolution evolution of galaxies 
with masses that go from dwarf galaxies {IQ^^Mq) to clusters of galaxies (lO^^M©). They fltted the final density 
profiles with typical parametrizations proposed by Navarro, Frenk and White [s^l (NFW): 



^(^) = / / ^2 ' (1) 



and by Moore [311 (M99): 



^^""^ (r/r,)i-5(^ + r/r,)i-5' 



finding that simulated data are well approximated by such profiles, that are, hence, "universal" , in the sense that 
the same analytical form successfully captures the shape of halos at different masses. However, the logarithmic slope 
(3{r) = d\np{r)/d\nr of the density profile decreases faster (more slowly) in the simulated data than does in the 
NFW (M99) profile at small radii. Moreover, iV-body simulations do not exhibit any indications that (3{r) converges 
to a central value /3o, as should happen for a NFW profile {(3q = —1.0) or for a M99 profile (/3o — —1.5). This can 
be due to the finite resolution of numerical simulations, which can be trusted down to the resolution radius rmm, 
usually taken to be around 0.5% of the virial radius, depending on the total number of particles in the simulation 
{Tmin ~ kpc, for Milky Way-sized halos). The structure of the inner region therefore is not clear yet, and the value 
of the central slope /?o can only be inferred by extrapolation. Since the region near rmm is where the deviations from 
Eqs. [1] and [5] are stronger, the extrapolation procedure can lead to significant errors. 

Aside from these uncertaintes, Navarro et al. [2^ used the density of their simulated halos to state that, given 
a generic power law profile whose central slope converge to a finite value, such value cannot be steeper than —1.5. 
Throughout this paper, we will consider this value as the lower limit for the slope of a DM nuclear profile in absence 
of BH. The profiles in Eqs. [T] and [2] are particular cases of a more general parametrization, the so-called {a,(3,j) 
profile [32: 

^('^)=^°(7)'[l + (r/ro)"](^-'r)/.' (3) 
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that reduces to p oc (p oc r~^) in the hmit of small (large) radii, while a characterizes the sharpness of the change 
in the logarithmic slope. The NFW is recovered for (a, (3, 7) = (1, 3, 1) and the M99 for (a, /3, 7) = (1, 3, 1.5). 

It has been proposed that, for a choice of parameters mimicking the presence of flat cores (e.g. if (a,/3,7) = 
(2, 3, 0.2)), such profiles provide good fits of both the rotation curves from direct observation of LSB galaxies and of 
the rotation curves from simulated nuclei, suggesting that core profiles should be preferred to Eqs. [T] and [2] since no 
contrast between A^-body techniques and observational data is present [s^. A possible explanation for the different 
results obtained by Navarro et al. [281 and Kravtsov et al. ^33*1 is the fact that the analysis by Navarro et al. does 
not take into account the large uncertainty in the determination of the luminosity distance of galaxies. 

Moreover, it was recently shown [s^l that the best fit to simulated data for high-resolution ACDM halos is obtained 
with profiles inspired from the so-called Sersic law [s^: 

ln(S]/Ee) = -6(Xi/"-l). (4) 

Such a relation provides the best description of the luminosity profiles of elliptical galaxies and the bulges of disk 
galaxies [3^: E is the projected density, X = R/Re, and R is the projected radius. The parameter n, called Sersic 
index, defines the shape of the profile, and 6 is a function of n, usually chosen so that the radius Re contains half of 
the luminosity of the galaxy. 

Eq. [Hcan be re- written as 



dlnE b f R 



d\nR n \ R, 



l/n 

(5) 



making explicit the power law behaviour of the logarithm slope. Parametrizing the spatial DM profile in a similar 
way, we obtain 

dlnpir) ^sf^V^" (6) 



d\nr \f-2 

or p{r) « exp(— ^r^/"), that is called Einasto profile^ in order to emphasize the difference from the Sersic law, since 
we are using now spatial, not projected, quantities. r_2 is the distance from the center where the slope is equal to 
—2. Eq. [S]was tested fitting the density of the DM halos simulated by Navarro et al. [2^, providing better results 
than a NFW profile, at least for dwarf and galaxy halos [l^l • The values of n (now called the Einasto index) , left as a 
free parameter in the fit, range from 4.33 to 7.44. A more detailed description of this model can be found in Merritt 
et al. [37| , where fits to simulated DM halos are studied and compared with alternative parametrizations. A number 
of recent studies [ssl . [39| have confirmed that simulated DM density profiles deviate slightly but systematically from 
the NFW form and are better approximated by Einasto's empirical law. 

If the Einasto relation was confirmed as a good parametrization of the inner region of DM halos, this would suggest 
that a scale-free relation like Eq. [SI describing both dark and luminous matter, is a characteristic feature for systems 
that form via gravitational clustering. 



III. PARTICLE DENSITY AROUND ALREADY-FORMED BLACK HOLES 

A. The Fokker-Planck equation and the Bahcall-Wolf solution 

A population of particles (both stars and DM particles) around a BH can be described by a distribution function 
/(x, V, i), whose evolution is governed by gravitational encounters among particles [40| . In the small-angle approxi- 
mation, such distribution function slowly diffuses in the phase space (x, v) towards a steady-state configuration. The 
time needed to achieve this equilibrium solution is defined as the relaxation time t^ei. For the stellar population, 
assuming that all stars have the same mass poj : 



0.34(7^ 

trel 



G-^prrii, In A 
0.95 • 10^° yrs 



(7) 



200 kms"V VlO^^^opc-V \Mq J \ 15 



where a is the velocity dispersion, p the stellar density and In A, known as the Coulomb logarithm, comes from 
imposing a physical upper cut-off in the distribution of impact parameters for stellar encounters. In A is usually 
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related to the mass of the central BH (M,) expressed in units of stellar masses [41[: 



In A - In (^) = In ] = ln(iV./2). (8) 



\ 2Gmi J \ 2toj, 

trei depends on the distance from the center of the galaxy but usually, as in Eqs. [7] and [51 it is computed at the 
influence radius r^, defined as the radius at which the gravitational potential due to the BH is equal to the kinetic 
energy: 



GM, _ „ / M, 

V lO^M© J V2OO km s" 



- ^ " 11 T7^ ) ( . -1 ) • (9) 



In the case a singular isothermal density profile [40| {p{r) = a/2TTGr^), M{r < r/,,) = 2A/,; and this can be 
used as an alternative definition of the influence radius. For the Milky Way (MW), r/j « 3 pc according to both 
definitions. There seems to be a clear trend of relaxation times with the mass of the central BH, where smaller 
objects (corresponding to fainter nuclei) are associated to smaller relaxation times. It can be seen from Fig. I (taken 
from Ref. 23]) using the empirical M, — a relation p^ : 



M, = 5.72 • IO^Mq f J- ] . (10) 

VlOkms-V 

Nuclei can then be classified in two different categories. Those nuclei with a relaxation time larger that the Hubble 
time, cannot have already achieved their relaxed equilibrium configuration, so that their distribution will reflect the 
process of nuclear formation. They are called collisionless nuclei, and they are characterized by a central region with 
a low density of stars, since near the BH a core is present with a slope < 0.2 23i, at least for those nuclei where 
the influence radius is resolved. The "mass deficit" (compared to what one expects from the Sersic law ^]) is up 
to 4 times the mass of the central BH. There are, then, galaxies, like the MW and M32 that are characterized by a 
relaxation time smaller than 10^" yrs (at resolved radii, e.g. the MW has trei — 3.5 x 10^ yrs at radius w 0.1 r/i). These 
coUisional nuclei have already reached their steady-state configuration. Usually they are faint nuclei (My < —20) 
and, opposite to cores of collisionless nuclei, the innermost region exceedes the Sersic law, establishing an inner power 
law profile with slope steeper than « 1.5 or a compact stellar nucleus [i^ . 

In this Section we will restrict ourselves to the case of coUisional nuclei where the relaxation time is smaller 
than the Hubble time and the nucleus has, today, a relaxed, steady-state equilibrium configuration for the stellar 
population. The diffuse evolution of an isotropic distribution / is described by the Fokker-Planck equation [ioL [4ll|. 
where gravitational collisions are taken into account and parametrized as: 



4.86 



2 ,r..9f BFe d 



dt dE dE 
with 



-^EE%-B,^ 



(11) 



(12) 



q(E) f dE'f{E') + [ dE'q{E')f{E' 

J -00 J E 

De{E) ^ -GAn'^G^m^lnA dE'p{E')f{E'), (13) 

Je 



q{E, t)^l r^^^ v^r'^dr = \ I [2{E ~ <l))]'^/\^dr, (14) 



3 7o 3 

while p{E) — —dq/dE is the volume of phase-space accessible to stars with energy E. 

The equilibrium solution cannot be a Maxwellian distribution, since it would imply an unphysical stellar density 
near the BH [3] , given that stars cannot be present at radii smaller than the tidal radius rt , inside which tidal forces 
tear stars apart. As a consequence, the distribution function is set to zero for r < rt- 

The physical steady-state solution was determined by Bahcall and Wolf [1^ , following a previous work of Peebles 
[4^. They proposed that the equilibrium configuration is a zero-flux solution, and obtained a distribution function 
with a power-law behaviour f{E) oc with a corresponding power-law density profile p{r) — p^r^'^/^. They also 
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Figure 1: Relaxation times measured at the SMBH influence radius in the ACS/Virgo sample of galaxies (see also Cote et al. 
(2004)), versus the central stellar velocity dispersion. Filled symbols are nuclei in which the influenced radius is resolved. The 
star is the MW. Figure taken from Merritt (2006). 



numerically solved a Fokker-Planck-like equation, obtaining a profile that can be very well described by the zero-flux 
solution, in the inner region (r < 0.2 r^, where the cusp actually forms), and a Keplerian rise in the velocity dispersion 
a oc r^^/^. See Fig. 2 taken from Ref. [1^. Their solution has been confirmed by iV-body simulations [4l|, in which 
the assumptions of isotropy and small-angle, characteristics of the Fokker-Planck-like formalism, had been relaxed. 
The validity of Eq. [TT] relies on the following assumptions: 

• stars are point-like masses, described by a distribution function that evolves due to gravitational interactions 
with the central BH and among themselves. In particular, this means that encounters in which stars collide 
with each other are neglected 

• the small-angle approximation: the gravitational potential belongs to a particular class of interactions for which 
the net force experienced by a test particle surrounded by a population of other bodies with which it interacts, 
is mainly due to cumulative, weak encounters with particles far away, rather than to strong but infrequent 
interactions with close particles. If we call po the impact parameter for the interaction that causes a deflection 
of 45° in the velocity of the test particle, all the close encounters with impact parameters p < po count only for 
4% of the total net deflection [40|. Working under the small-angle approximation means that, in the study of 
the evolution of a test particle with a velocity v embedded in a larger particle population, we are considering 
only the encounters with those bodies far away that will produce small deflections Av compared to the initial 
velocity v 

• the distribution function does not depends on the angular momentum and, since the gravitational potential is 
function only of the radial coordinate, the stellar distribution is isotropic 

• the mass of the central BH is much larger than the mass of a star and does not change with time, so that inside 
the influence radius, the gravitational potential is constant and Keplerian and is due only to the BH itself. 
Under this assumption we can (as in Eq. [TT]) neglect the term proportional to Of /dE that would appear in a 
more complete form of the Fokker-Planck equation. Moreover the BH mass has to be much smaller than the 
total stellar mass near the BH itself. Requiring this particular mass hierarchy [i^, leads us to a consequent 
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Figure 2: Evolution of the stellar distribution around a SMBH due to energy exchange between stars. These curves were 
computed from the isotropic, orbit-average Fokker-Planck equation with boundary condition / = at log \E\ — 6. Left panel: 
phase-space density /; right panel: configuration-space density p. The initial distribution (shown in bold) had p cx r""'^ near 
the SMBH; thin curves show / and p at times of (0.2, 0.4, 0.6, 1.0) in units of the relaxation time at the SMBH's initial 
influence radius Vh- Dashed line show the "zero-flux" solution / oc and p oc r~''^^. The steady-state density is well 

approximated by the zero- flux solution at r < 0.2 ru. Figure taken from Merritt (2006). 



timescale hierarchy, i.e., the assumption m^, <C M, <^ Mi,{r < rh) implies that the crossing time t^r (the time 
needed for a star to cross the nucleus) is much shorter than the relaxation time, so that in a crossing time stars 
do not experience any changes in the physical proprieties of the system. 

All the assumptions listed above can also be satisfied in the case of a distribution of DM particles only (with a 
common mass m^), so that a Fokker-Planck formalism is appropriate also for DM. We expect the existence of a relaxed 
solution also in this case, but the relaxation timescale for DM is enormously larger than for stars and DM particles 
will never reach their steady-state solution: they can be effectively considered as collisionless objects, pratically not 
sensitive to the gravitational self-interactions that drive the dynamical evolution described by the Fokker-Planck 
equation. 

The requirement of a common mass for particles (stars or DM) is not included in the list above, because in the more 
realistic case of a nucleus with particles of different mass, Eq. [TT] can be modified in order to describe a multi-mass 
case. For a two-component nucleus, made of stars (with a common mass m^) and DM particles with <C to* the 
Fokker-Planck equations will be the following: 



dt 



d_ 

dE 



-rrii^DEg* - Dee 



dE 



(15) 



Att'^p{E) 



dgx ^ d 
dt dE 



—Dee 



dgx 

dE 



with 



gi{E,t)= I fi{E,t,m)m dm, 
Jo 



hdE,t) 



fi{E,t,m)inn? dm, 



with i G x], and the diffusion coefficents can be written as 



Dee{E) = 647r*G''m*lnA 



q{E) r dE'g,{E')+ f dE'q[E')g,{E') 

J-oo J E 



(16) 

(17) 
(18) 

(19) 
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r 

De{E) ^ -GAn'^G^lnA dE'p{E')g^{E'). (20) 

JE 

fi, is the stellar distribution function, whose evolution (Eq. [T5)) is governed by star-star interactions, and is the 
distribution function for DM and in Eq. [16] only DM-star encounters are considered due to the coUisionless nature of 
DM. 

The final steady-state solutions will have the usual —7/4 slope for stars (in fact Eq. [15] is not different from Eq. 
[TT|) and a milder —3/2 slope for DM, that will be established in the same timescale trei (Eq. [7]). The steeping of the 
initial profile, due to the presence of a BH, leads to an increase of the DM annihilation rate. This is why we refer to 
BHs as DM Annihilation Boosters. Although in this case the effect is not dramatic, since a r~^'^ profile can be found 
even for models without BHs (see Sec. [IT]), and it is more likely that the overdensity will be reduced in a couple of 
relaxation times [l^l (see Sec. IIV C|) . we will see cases in the next Sections where BHs can provide huge boost factors. 

The presence of a stellar cusp has been experimentally confirmed for the MW ^il, , through the detection of 
a profile with a slope equal to —1.4 < 7 < —1.3 in the inner region (r < 0.38 pc) and to —2 (isothermal profile) 
in the outer region (see Fig. 3). Our Galaxy is a collisional nucleus, since the relaxation time is shorter than the 
age of the Universe (3.5 x 10^ yrs at ~ 0.1 r^), so it was suggested to interpret its cusp as the Bahcall-Wolf solution 
to the presence of a SMBH with a mass « 3.7 x 10® M© (r/i is « 3 pc so that the cusp starts more or less where 
r = 0.1 rfi), hypothesis supported also by the luminosity of the X-ray source Sgr A*. Anyway, it is more likely that the 
MW experienced a merger between a redshift z = 2 [50| and today, so the cusp will be the result of an "overdensity 
regeneration" (see Sec. IIVC|) . The actual, detected profile [5l| is consistent with a cusp regenerated after a merger 
occured at a time > 8 Gyr in the past. 

The detected inner slope of —1.4 is not exactly what the Bahcall-Wolf solution predicts (7 = —1.75). However, 
the two results are considered consistent with each other since the steeper value is derived under the simplyfing 
assumption of a population of stars with identical mass, and if the more realistic multi-mass formalism is introduced, 
the slope will become shallower, moving towards the —1.4 value. The same can be said if, as it was argued [5l|, the 
time required to reach a steady-state solution at the Galactic center is > lO^^yr. 
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Figure 3: Mass density profiles near the center of the MW and M32. Dashed line is p(r) oc r Both galaxies contain SMBHs 
with masses ~ 3 • 10® Af© and with infiuence radii rn ~ 3 pc. Figure taken from Merritt (2006). 
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B. Loss-cone dynamics and BH binaries 



In this Section, we focus on two mechanisms that can reduce, and even completely destroy, coUisional cusps. First 
of all, we review the basics of loss-cone dynamics, i.e. the effects related to the presence of a BH tidal radius or of the 
relativistic event horizon. This is only partially relevant for DM since it will affect mainly the stellar population, but 
it is relevant in this contest since we have just seen (Eq. [T6|) how stars and DM evolve together, so that a decrease 
in the stellar cusp due to loss-cone \23\ will influence the DM density, too. Moreover we want to emphasize that any 
enhancement in the density of DM and stars has to survive to a collection of dynamical mechanisms that can damp 
them. 

After discussing the loss-cone, wc will focus on the possibility that, during the merging evolution of a nucleus, a 
BH binary forms. In this case, both stellar and DM distributions will be directly affected leading not only to the 
reduction of cusps but also to the possible destruction of spikes in models with adiabatic growth [sO] (see Sec. lIVBp . 

The loss-cone of a BH is the set of eccentric orbits populated by stars that are going to intersect the tidal radius. 
Such orbits are depleted in a crossing time icr, since stars are eaten up by the BH: tidal forces inside the tidal radius 
tear a star apart and these events are accompanied by flares, emissions of light peaked in the X or UV band with a 
luminosity of ~ lO^^'erg s~^. A handful of these flares have been observed: they have the expected signature and the 
number of detections is roughly consistent with theoretical estimates of the consumption rate N jssllM. [ssl . [56j . 

Once depleted, loss-conc orbits can continue to contribute to the consumption rate only if they are somehow refilled 
with particles. Energy diffusion provides a mechanism for such a repopulation. In fact, A^-body simulations confirmed 
[23I I that the zero-flux Bahcall-Wolf solution is established only approximately and that (for stars) there is a residual 

This value is too low compared to the expected N whose main contribution comes, instead, from angular momentum 
diffusion: the "classic loss-cone theory" (4^,113,1111 apphes to globular clusters (populated with a central BH), whose 
relaxation time is so low that they are well-relaxed and old objects. The distribution of stars near the BH is therefore 
assumed to be on a steady-state, for which a Fokker-Planck-like formalism is appropriate. Resulting estimates for the 
consumption flux can be introduced in the original Fokker-Planck equation (Eq. to study how particle density is 
affected by the presence of a tidal radius: 

4nME)^--^+PiciE,t). (21) 

When extended to the study of galactic nuclei, the steady-state approximation may fail, at least for coUisionless 
nuclei, and the stellar profile near the tidal radius is, in general, different from the Bahcall-Wolf one. For example, 
galaxies are only approximately spherical: their shape is more likely to be triaxial and there is the possibility that 
they are governed by centrophillic orbits, i.e. orbits that pass arbitrarily close to the BH. In the case that these 
chaotic orbits survive until late stages in the galactic evolution, A^ would increase, since more particles would fall into 
the tidal sphere. 

Alternatively, the present galaxy can be the result of cumulative mergers of less massive mini-galaxies, each of them 
hosting a mini-BH: the formation of a BH binary would decrease the consumption rate, since all stars with angular 
momentum L < Ltin = {2GMi2ah) would be ejected [1§| (M12 is the total mass of the binary and Uh the major 
semi-axis when the system becomes "hard" ) , preventing loss-cone repopulation and leading to lower rates N. 

Finally, in real galaxies, the diffuse mechanism of orbits refillment will cause the nucleus to expand [2^ [goI. [giI. [6^ . 
since the density is reduced when particles are eaten up and those particles which fall into the loss-cone transfer energy 
to the remaining nucleus with the same effect of a heating process. The expansion is visible in one single relaxation 
time, the "decay" goes on at a constant velocity and the density can be written as p{r,t) = pc{t)p*(r), where p*{r) 
is the initial profile, while pdt) oc t~^. As a consequence, present-day collisional nuclei could have been denser in the 
past. 

It has been suggested that also the presence of a BH binary can effectively reduce the Bahcall-Wolf cusp. The 
growth of a galaxy is thought to pass through the agglomeration of smaller galaxies and protogalactic fragments. If 
more than one of these subhalos contain a BH, the two objects will form a binary system whose dynamics can strongly 
affect stars and DM. This scenario has received great attention since mergers and the ultimate coalescence of the BH 
binary are ideal targets for the detection of gravitational waves (63j . Evidences for the presence of such binaries can 
be found in Ref. [6^ and are based on the detection of multiple active nuclei in the same galaxy [65| . 

Consider a compact object with mass M2 moving, with its nucleus, around a BH with mass Mi, being q = M2/M1 < 
1 the mass ratio and M12 the total mass. The evolution of the binary can be described by three different phases 
[H, [13|: first, the smaller BH decays due to the dynamical friction with stars of the other nucleus, and the separation 
i?i2 between the two objects drops down. When the influence radius of the more massive BH is reached, the two 
objects can be considered as a bound object and the first phase comes to an end. The infall time scale [11] suggests 
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that binaries are not so uncommon since, for reasonable values of g (q « 10~^ and Mi « 10^ Mq), they are able to 
form before the Hubble time. 

The second phase is characterized by a quick "shrinking" of the binary, until it becomes hard, i.e. the binding 
energy equals the kinetic energy, or equivalently the major semi-axis reaches 52] — q/{l + g)^ • r/i/4. The third 
phase, when a < a^, is the least known: a binary in a fixed background begins to harden at a constant hardening rate 
s = d{l/a)/dt, but physical binaries has already ejected almost all stars on intersecting orbits and the rate suddendly 
drops. These orbits need to be repopulated, usually by energy diffusion, but this effect is more likely to be only 
subdominant, at least in those bright galaxies where the scouring of BH binaries has been detected, characterized by 
a relaxation time higher than the Hubble time. In numerical simulations with finite N, gravitational encounters will 
unphysically continue to supply particles to the binary at rates roughly proportional to N, so experimentally it is 
more usefuU to define the semi-major axis a^taz; where the hardening rate goes to zero [s^ . From A'^-body simulations 
it results 



dstall 



r 



is a second influence radius, defined as the radius where the total mass of particles within after the binary 
has stalled is equal to twice M12 . This values for astaii is a couple of orders of magnitude higher than the distance 
where the binary coalescences. This is known as the "final parsec problem" , since evidence is strong that BHs binaries 
do efficiently coalescence [H, at last. Many solutions have been proposed [G^I, e.g. that, as for loss-cone orbits, 
the presence of centrophillic orbits in realistic triaxial galaxies can affect the above considerations, so that the binary 
continues to shrink even to a ah- 

Assuming that the binary does stall at astaii, it will have transferred an energy 

Ai? « + ^MlMl « ^ 2(Mi + M,)a^ « 2(Mi + M,)a\ (23) 

2rh 2ah 2 

to the particles in the nucleus. This relation has been used to explain the mass deficit in the core of brightest 
galaxies since such an energy release will let stars leave the central core, with a total displacement of mass Mdef and 
0.4 < Mdef/Mi2 < 0.6 ^ for 0.05 <q< 0.5. 

Strictly speaking, the observed mass deficits M^ef reach values that are even four times larger than the mass of the 
binary M12 (that, if the coalescence occurs, is also the mass of the final BH). We can account for values as large as 
Mdef /M12 ^ 2 if the nucleus experiences more than one merger, with more than one binary forming. The total mass 
displaced will be simply the sum of each M^ef during each single merger. For even larger values, other mechanisms 
have to be evoked, e.g. the possibility that a third BH arrives when the first two have not coalescenced yet. In such 
a situation, one of the BH usually leaves the nucleus {gravitational slingshot effect) leading to higher values for Mdef- 
Similarly one of the SMBH of a BH binary can be expelled with a high velocity, due to the so-called gravitational-wave 
rocket effect [6§| . 



IV. ADIABATIC GROWTH OF BLACK HOLES 



A. Adiabatic growth of Black Holes 



In this Section we relax the assumption of time-independent gravitational potential. As we will see, this can, in 
some cases, lead to large DM overdensities. In particular, the adiabatic growth of BHs can produce the steepest DM 
profiles discussed in literature. 

The seed BH grows in an already-formed nucleus with a stable (stellar or DM) configuration. The condition 
of adiabaticity guarantees that the growth timescale is large compared to the crossing time, but smaller than the 
relaxation time (for the stellar distribution). As a consequence, nuclei where a BH have grown adiabatically have 
not yet reached a stable, relaxed stellar configuration. Another consequence (69j is that the angular momentum and 
the radial action (i.e. = f Vrdr, where Vr is the radial velocity and the integral is over one closed orbit) are 
conserved. The hierarchy between the BH accretion timescale and the nuclear crossing time, that lies at the core of 
the adiabatic assumption, is reasonable at least for BHs with masses M, < lO^^M©, as can be checked adopting the 
shortest timescale for the BH growth, i.e. the Salpeter time ts = M,/MEdd (where MEdd is the Eddington accretion 
rate), and comparing it with the crossing time at the influence radius tcr oc M,a/rfi, using the M, — a relation (See 
Eq. 

The first study on the impact of adiabatic growth [70| on stars, analyzed the case of a non-singular isothermal 
stellar profile, and predicted an overdensity extending to the same size of the initial core, with a slope equal to —3/2. 
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In the case of a DM halo, such overdensity has been called spike [7l|, to distinguish it from the aforementioned DM 
cusps. 

A numerical algorithm that mimic adiabatic growth was also developed, in order to confirm the creation of the 
overdensity |72| . The method is very flexible and, in fact, it was applied to initial models other than the isothermal 
distribution |72 |. Two classes can be identified: the first includes all those profiles called "analytic cores" , characterized 
by a density that can be expanded in a power law series near the BH (p(r) w po + l/2/3oV^ + ...), while the second 
describes the so-called 7 models that exhibit a power law density profile in the inner region: p{r) oc r~^. 

As benchmark cases, the 7-models with < 7 < 2 and the isothermal model (as an example of analytic profiles) 
are considered here, and results are presented in Tab. I and in Fig. 4 [t^ (only for the isothermal and for 7 = 1). 

The spike radius r^p, i.e. the distance where the slope changes due to the presence of the BH, depends on the BH 
mass and it is related to the influence radius as r^p « 0.2 |74| . Inside such radius, the spike has a slope 73^, that 
depends on the initial 7. In the case of a model with analytic core the final slope is —3/2 70], while for the 7 models 
an analytic relation holds [iej [t^-ItsI: 

Is, = (24) 

4 — 7 

Such relation is valid under the following assumptions (satisfied by all the models in Quinlan et al. [72i] ) : 

• the distribution function is isotropic; 

• the gravitational potential can be written as r^^^ in the r — > limit; 



/ diverges as [E — $(0)] " in the limit E ^{0) (this last requirement is what makes a model with analytic 
core different from a 7 model). 



Model 


7 n 7sp C 


isothermal 
7 model (7 = 0) 
7 model (7 = 1) 
7 model (7 — 3/2) 
7 model (7 = 2) 


3/2 9/4 

1 2 9/4 

1 5/2 7/3 7/3 
3/2 9/2 12/5 12/5 

2 5/2 5/2 



Table L Different quantities computed from the adiabatic growth of the initial models proposed by Quinlan, Sigurdsson and 
Hernquist (1995). 7 and 7sp are the initial and final slope in the density profile for the region closer to the BH. n indicates 
how the distribution function diverges as _E ^ '1*(0) and C is the slope of the final density profile if it were made of particles 
on circular orbits. The value for n in the 7 = 2-model is absent since the equation used to derive n is not valid for 7 = 2, but 
for 7 — > 2 the final profile has ysp 5/2. 

Comparing the first two lines in Tab. I, it can be seen that, even if both models start with a costant core, they 
develop very different final spikes, due to the different behaviour of the distribution functions in the E ^'(0) limit, 
suggesting that the formation of a strong spike is not a consequence of a singularity in the density profile but in the 
distribution function, and in particular in the way cold orbits (populated by stars with a low velocity) are arranged 
(see Sec. IIVBI) . We will not consider initial configurations with 7 larger than 1.5 (see Sec. [TTI, so the steepest spike 
has a slope of —12/5 (when 7 = 1.5). 

The velocity dispersion reacts to the BH growth in a similar way for both classes of initial models: in fact, a 
Keplerian rise appears, with a slope of —1/2. On the contrary, the anisotropy is substantially different: analytic 
models exhibit a mild tangential anistropy at an intermediate distance from the BH but remains isotropic at the 
center; the more massive the BH is, the higher the anistropy. While for 7 models, orbits are tangentially-biased in the 
central region and the anisotropic area increases with more massive BHs. If the hypothesis of an isotropic distributions 
is relaxed, a nucleus made interely by circular orbits evolves to a profile with final slope equal to C, shown in Tab. 
I. As one can see, for 7 models 7^^ = C, while the circular isothermal model exhibits a much steeper slope, although 
with a value not higher than for the 7 models. Such consideration suggests that results from adiabatic growth are 
not very sensible to possible violations of isotropy in the initial configuration. 

Even if spikes are the steepest known overdensities, they cannot be considered as signatures of BH growing adia- 
batically (since a simple singular isothermal profile that nothing has to do with adiabatic growth is steeper than half 
of the models in Tab. I) and neither hints of the presence of a central BH (since a —1.5 slope, as in the first line of 
Tab. I, can be found for halos without BH). On the contrary, a rise in the a profile is quite a robust indication of the 
presence of a central object. 
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log(E) log(R) 

Figure 4: Adiabatic growth in the isothermal {left panel) and for a 7 = 1 {right panel) model for a value of Af, of 0.001, 0.01, 
0.1 and 1.0 (mass units explained in the Reference), with the mass encreasing from bottom to top in the top two panels and 
from top to bottom in the last panel. The density is shown in the top frame, the averaged radial velocity in the middle and 
the anisotropy parameter f3 in the last frame. /3 — 1— {vt)/{2vr) where vt {vr) is the tangential (radial) component of velocity. 
Figure taken from MacMillan, Henriksen (2002). 



B. Destruction of spikes 

The formation of spikes described in the previous Section leads to the largest annihilation boost factors: the 
mechanism of adiabatic growth can, in fact, produce inner slopes as steep as —2.25 (See Tab. I), for profiles that will 
be characterized by a large annihilation flux (which is proportional to the integral of the DM density squared) , with 
interesting consequences for Indirect DM searches. 

However, it has been argued [t^ [IOjIzEIi that the formation of spikes requires ad hoc initial conditions for the DM 
halo. Moreover, even if spikes do form, then, dynamical effects can reduce or even destroy them, as we considered 
before for coUisional cusps. A spike can form even from an initial density profile that does not diverge [t^, but in 
order to produce a significant overdensity, the distribution function of cold orbits has to diverge in the E ^ ^{0) 
limit. In fact, cold orbits are those which provide the particles that will form the spike, since they are the most 
affected by the presence of the central BH. 

But these cold orbits are more likely to be depleted due to the interactions of stars with molecular clouds or globular 
clusters or other bodies that can pass through a galactic nucleus. Moreover, the evolution of a galactic nucleus is 
thought to be characterized by the cumulative mergers of sub-nuclei and even a single merger event can have dramatic 
consequences on the distribution of cold orbits. In other words, one can enumerate a collection of effects that effectively 
heat up the particles near the BH, so that thay can leave the central region obstacling the formation of the spike. 

Also in the case that the spike is formed, it is unlikely that it will survive to the evolution of the nucleus and, 
in particular to the presence of dynamical mechanisms that would provide an additional heating source to particles 
on cold orbits, with the result of highly reduce of even destroy the enhancement. Numerical simulations have been 
performed in order to quantify these effects: for example the possibility that the BH forms slightly off the center of 
the nucleus was described by Ullio, Zhao and Kamionkowski [73 ■ The BH would slowly spiral in, towards the center 
fr^ . and then adiabatically grow to the final value. But if the initial value for the BH mass is too low, the spiraling 
would take too long to finally reach the center, while, if the BH is too massive, its scouring effect on the DM particles 
would flatten the central density, to values that can be even lower than the initial profile (see Fig. 5). 

Moreover, gravitational interactions of DM particles with baryons in stars modify the evolution of DM in the spike, 
reducing the enhancement, in the same way that stars heat the DM particles in a coUisional cus p c ausing its damping 
(see Sec. IIV CI) (47j (75| . Simulations on the effects of galactic mergers can also be found in Rcf. [soj . Other objections 
have been put forward, suggesting that spikes can form only as results of a series of accidents and, therefore, are not 
expected to be common in the local Universe [75i] . 
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Figure 5: Modification of a NFW DM density profile due to tfie off-center formation of a BH seed of mass M§u, its spiraling 
towards tfie center of tlie DM distribution, and its adiabatic growth to the present-day mass of the BH Mbh at the Galactic 
center. Several different values of the BH seed mass are plotted, pcore is the maximum WIMP density above which WIMPs 
are depleted by pair annihilations. Figure taken from UUio, Zhao and Kamionkowski (2001). 



C. Regeneration of cusps 

In the previous Sections we focused on processes that can destroy overdensities (whether cusps or spikes). However, 
gravitational interactions among particles (during the evolution of the nucleus) can partially regenerate such structures. 
In a realistic nucleus populated by stars and DM, star-star collisions and star-DM collisions (given enough time) drive 
the evolution towards a steady-state —7/4 profile for the baryonic component and a —3/2 profile for the DM component 
[131 (see Section IIII Ap . This is true also if the nucleus is the result of an early evolution phase in which previous 
enhancements were destroyed. In other words, consider a nucleus with a short trei and a steep DM and stellar profile 
(due to coUisional relaxation or adiabatic growth of the central BH). If a merger occurs and a BH binary forms, the 
displaced mass will reduce or even destroy both overdensities. But, due to the short relaxation time, it can happen that 
the nucleus has enough time to reconstruct, from the core profile after the BHs coalescence, the collisional solutions. 
The new DM cusp is called CREST (CoUisionally REgenerated STructure). This idea can be checked analitically, 
applying the two-body Fokker-Plack formalism to a core profile describing a nucleus after the scouring of a BH binary, 
but also with the help of numerical routines [IJ. Results are summarized in Fig. 6: DM CRESTs are not as steep as 
spikes, but they have the advantage to form from very general initial profiles, given that the (stellar) relaxation time 
is short enough. They need a timescale of roughly trei{0-2 r^) to form, but then DM particles in the CREST continue 
to be heated by gravitational interactions with stars and the —3/2 solution, therefore, decays in a self-similar way 

Px(r,t) «p^,oWG(i/irei), (25) 

with dG/dt < 0, so that after 4.5 trei{Th) the reduction is of a factor 1/e^. 

The balance between the requirement that the relaxation time is short enough to let the CREST form but not too 
short to make the CREST not to decay too much, leaves us with a rather narrow window of galaxies where CRESTs 
can be present: one can detect them in galaxies with a luminosity 3 • 10*^0 < L < 3 • lO^L©. The MW is inside 
this range and, in fact, many proposed to interpret the cusp detected for our Galaxy [49| as a reconstructed structure 
after a merger occured ^ 8 Myr ago [5l| . 
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Figure 6: Solutions of the Fokker-Planck equation that describe the joint evolution of stars and dark matter around a BH due 
to star-star and star-DM gravitational encounters. Length unit is rs, while the density is in units of the initial value at rh- 
Curves show the stellar (left) and DM (right) density profiles at time 0, 0.2, 0.4, 0.6, 0.8, 1.0 in units of the initial relaxation 
time at rh- Dashed lines are the steady-state solutions. Figure taken from Merritt, Harfst and Bertone (2006). 



V. INTERMEDIATE MASS BLACK HOLES 



A. Formation scenarios 



The steep DM slopes produced by the adiabatic growth of BHs make these objects extremely interesting since they 
may effectively act as DM Annihilation Boosters. This circumstance encouraged many authors to look for possible 
ways to evade the dynamical effects causing spikes to damp. 

Two possibilities have been proposed: the first is to focus on IMBHs instead than of SMBHs. IMBHs can be present 
within substructures of DM halos [l^, [2l| and their evolution is such that the objections raised for spikes around 
SMBHs do not apply [22] . The second is considering the contribution of spikes and mini-spikes (DM overdensity around 
IMBHs) to the diffuse Extragalactic Gamma ray Background (EGB), integrating the signal from high redshift, i.e. 
when spikes were already formed but destruction mechanisms were not yet effective [tI, [t^ [S^I ■ 

We focus on IMBHs first. Their masses range from w 20 to IO^Mq and adiabatic growth leads to the formation 
of "mini-spikes" in the same way as with SMBHs. We don't have direct observational evidence for the existence of 
IMBHs, but some hints come, e.g., from Ultra Luminous X-ray sources (ULXs) [8l|, sources that emit in the X band 
with a luminosity higher than lO'^^erg s''^, and, hence, not compatible with the interpretation as BHs accreting at 
the Eddington limit. But, due to their positions in the host galaxy, they cannot be explained in terms of AGNs 
either. The hypothesis of a BH with a mass higher than 15 — 20Mq and less massive than a SMBH seems to be a fair 
explanation, instead. 

Many authors also proposed that globular clusters can host IMBHs, and a possible confirmation of this hypothesis 
comes from the fact that the mass scale for an IMBH and the value of the stellar velocity dispersion measured in 
globular clusters fall exactly at the extrapolation at lower values of the M, — cr relation valid for SMBHs |20j| . From 
a theoretical point of view, IMBHs can also help to explain the formation of SMBHs: the Sloan Digital Survey [s^ 
[S^ . [s^ l has detected quasars up to redshift z « 6 suggesting that SMBHs were already present when the Universe was 
~ 1 Gyr old. One of the most natural way to understand this is that SMBHs grew, through a phase of fast accretion 
and mergers, starting from already massive seeds. In fact, a generic prediction of scenarios that seek to explain the 
properties of the observed SMBH population, is that a large number of "wandering" IMBHs exist in DM halos [85l.[86j 

Despite their theoretical interest, it is difficult to obtain conclusive evidence for the existence of IMBHs. A viable 
detection strategy could be the search for gravitational waves produced in the mergers of the IMBHs population 
[6^ . [stI [qoI . [91| , with space-based interferometers such as the L arg e Interferometric Space Antenna LISA [92| . 

Two formation scenarios are discussed here, follwing Ref. [22|. In the first (scenario A), IMBHs form from the 
gravitational collapse of Population HI stars, which are usually heavier than local stars, since they grow in an 
environment with very low metallicity, for which metal line cooling can be neglected. As a consequence, the Jeans 
mass (that scales with the temperature as T'^/^) is higher, allowing the formation of more massive structures. Such 
stars are characterized by very low metallicity, too, meaning that they will lose little of their mass due to winds and 
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weak pulsations. Population III stars with masses above larger than 25OM0 would be able to collapse directly to BHs 
without any explosion |i2Qi] . 

In the second scenario (scenario B), the formation starts at high redshift (z « 15) from a gas cloud that is massive 
enough to cool down forming a protogalactic pressure-supported disc at the center of the cloud, made mainly by 
stars that, according to the distribution of angular momentum in the cloud, are in the low momentum tail of the 
distribution. The dynamics of the disc is governed by an effective viscosity that transfers matter at the center and 
angular momentum outward. This flow will go on until the first stars start to fragment in the outer region of the disc: 
the so-formed central object will undergo gravitational collapse, forming the final BH, with a mass scale of » 10^ Mq 
(corresponding to an initial cloud with virial mass Myir ~ lO^M©) or [22l | 



M. = 3.8 X 10%/e f^) r^V'' (^) X f^V^'^ ( ^) ^ (^6) 

where / is the fraction of the total baryonic mass in the gas cloud that has cooled into the disc, k is that fraction 
of the baryonic mass of the disc that forms the final BH, M^r is the halo virial mass, z/ is the redshift when the 
formation starts from the cloud and t the timescale for the evolution of the first generation of stars. 

Regardless of the particular formation scenario, a population of IMBHs is predicted, each of them surrounded by 
a mini-DM halo and a mini-spike. These mini-halos evolve, merge with each other and form the actual big halos of 
galaxies. The hosted IMBHs will merge too, possibly contributing to the formation of the central SMBH. However a 
fraction of initial mini-halos never experiences any mergers and the associated DM overdensity can survive till today. 

Moreover, at least for scenario B, IMBHs form exactly at the center of the baryonic distribution in the mini-halos, 
forbidding BH off-center formation. Strictly speaking, in order to avoid off-center formation, one should require that 
the pristine BH forms at the center of the DM distribution and this does not necessarily coincide with the center of 
baryons, since stars, been collisional, can experience a different evolution than the coUissionless DM, resulting in a net 
displacement between the two distributions. But mini-halos are supposed to have a very low baryonic content, with 
no violent interactions able to drive the two distributions away one from the other: the fact that mini-halos are made 
most entirely by DM solves the possible off-center formation and, at the same time, the problem related to stars-DM 
interactions. Scenario B is thus able to circumvent all the mechanisms for spike destruction described above, while 
in scenario A they are reduced but still efficient: in the following Sections, we are going to consider scenario A as a 
conservative model for IMBHs, instead of the more optimistic scenario B. 

For the MW, an IMBHs population of 1027 ± 84 objects is predicted by the numerical simulations in Ref. [2^ . 
with a mass of approximatively 10^ Mq (scenario A), while for the scenario B the MW hosts just 101 ± 22 with a 
distribution in mass centered on 10^ Mq and log-normally distributed (cr, = 0.9), as it can be seen in the left panel of 
Fig. 7. 



B. Gamma-rays from DM annihilations around IMBHs 

The prospects for indirect detection of DM from mini-spikes around IMBHs populating the MW halo has been 
studied by Bertone, Zentner and Silk [l^. They considered 200 different statistical realizations of the MW, obtaining 
an IMBH catalogue where each object is surrounded by a mini-DM halo with an inner slope of —7/3, resulting from 
adiabatic growth of an initial NFW profile (Eq. [1]). 

The spike extends from a cut radius Vcut to the outer spike radius r^p w 0.2 r/j. The cut radius depends on the 
mass and annihilation cross section of the DM candidate, being defined as the radius where the DM density reaches 
an upper limit due to DM annihilations [22]. 

Eq. [27] was used to compute the annihilation flux from each mini-spike: 

^=rdE'^ (27) 
la.J_diV, p 

dN~t / Gv \ / TOv f d 



dE Vl0-26cm3s-iMl00 GeV/ Vkpc 
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Figure 7: Left panel: Mass function of unmerged IMBHs in scenario B, for a MW halo ai z = 0, averaged over 200 Monte Carlo 
realizations of a MW-like halo with virial mass M^tr — W^^'^ Mq. Figure taken from Bertone, Zentner and Silk (2005). 
Right panel: Differential photon spectrum per annihilation. Different parametrizations and annihilation channels are shown. 
The solid line is an analytic fit to the 66 channel, as obtained by Fornengo, Fieri and Scopel (2004) with = 1 TeV. The 
dashed line is relative to the same annihilation channel, but with a different parametrization of the Fragmentation Functions 
obtained by Fornasa, Taoso and Bertone (2007). The dotted line is relative to UED DM and includes FSR from annihilation 
to charged leptons. Figure taken from Fornasa, Taoso and Bertone (2007). 



where $0 = 9 • 10^^"cni^^s^^, is the mass of the DM particle and av its thermally-averaged annihilation cross 
section, d is the distance of the IMBH from the detector and dN^/dE is the number of photons with energy from E 
to E + dE produced by the annihilation. The differential flux is then integrated in energy from a lower threshold Eth 
that depends on the experiment considered, up to m^. 

The differential spectrum dN^/dE can be written as a sum over all possible annihilation channels of the number 
of photons dN^/dE produced from the particular channel xx ^ weighted by the corresponding branching ratio 
S" = BR(xx^aa): 



dE ^ dE 



B^-^. (28) 



a 



If DM particles annihilate directly into photons, the spectrum will be a line at an energy equal to the DM particle 
mass, but this contribution is usually suppressed at least for a neutralino-like DM candidate (a branching ratio to 
primary photons of « 10"'^ is used in Ref. 79]). Otherwise the spectrum will be a continuum, and we can distinguish 
between two classes of spectra: in the first one photons are produced from the decay of neutral pions (tt" — > 77) 
formed through the fragmentation of annihilations products such as heavy quarks, gauge or Higgs bosons. Many 
parametrizations are available for these spectra dN^/dE [9^ [§3, IsH], that represent a suitable approximation for 
neutralino-like candidates. In the right panel of Fig. 7, two possible examples for the fragmentations of the neutral 
pion 7r° are compared [qI, [9^1 , in the easy-to-consider case of DM annihilating only to a pair of b quarks . 

The second possibility is inspired from candidates arising in theories with Universal Extra-Dimensions, where 
the DM can annihilate to charged leptons with a large branching ratio (for neutralinos, such channel is severely 
suppressed). Such final states, then, can produce photons through the so-called Final State Radiation (FSR). A 
typical FSR spectrum is plotted in Fig. 7 with a 40% of DM annihilating to fight fermions, 20% to rs and the 
remaining in b quarks. The photon production through light fermions is computed analytically [9^ while the other 
two channels are described as in the NPD spectrum. 

Differences between the two classes of annihilation spectra become more pronounced above x « 0.6, since at low 
energy the pion decay dominates. The FSR spectrum is characterised by a flatter spectrum (the slope in approximately 
— 1 instead of —1.5 for NPD decay) and by a more abrupt cut-off at an energy equal to the DM particle mass. 

The prospects for detecting IMBHs in the MW are summarized in Fig. 8 where the number of point-like sources 
(associated to IMBHs) with an annihilation flux higher than $, is plotted as a function of <&. Compared with the 
sensitivities of GLAST [93] (a Space-Based telescope scheduled to be launched at the beginning of 2008) and EGRET 
[gst (a Space-Based gamma-ray detector, whose data are available on Ref. [9^) for a 5a detection, the most optimistic 
configuration (m^^ = 100 GeV and av = 3 ■ lO^^^cm^^s^^) leads to almost 100 (80) detectable sources by GLAST 
(EGRET) in 1 year. 
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Figure 8: IMBHs integrated luminosity function, i.e. number of BHs producing a gamma-ray fiux larger than a given flux, 
as a function of the flux, for scenario B and a NPD spectrum. The upper (lower) line corresponds to — 100 GeV, 
av = 3 ■ 10~^^cm^s~^ (m^ = 1 TeV, av = 10~^^cm^s~^). For each curve we also show the la scatter among different 
realization of the MW DM halo. The figure can be interpreted as the number of IMBHs that can be detected from experiments 
with a point-source sensitivity $ (above 1 GeV), as a function of We show for comparison the 5cr point-source sensitivity 
above 1 GeV of EGRET and GLAST (1 year). Figure taken from Bertone, Zentner and Silk (2005). 



The prospects for Indirect DM detection usually strongly depend on the particle physics parameters. For IMBHs, 
however, the av/m^ dependence of the flux in Eq. [27] is modified by the implicit dependence of the cut radius on the 

mass and the annihilation cross section av, so that at the end d^/dE cx ((tu)^/''tox^^^. 

In order to discriminate IMBHs from other extragalactic sources, one can perform a maximum likelihood analysis 
of the annihilation spectrum or look for a class of sources with an identical cut-off, at the mass of the DM particle 
[2^ . Alternatively, one may look at the Andromeda galaxy M31 fosl], which is located » 780 kpc away from us, and 
where a population of 65 ± 15 objects (in scenario B) is predicted. A fraction of these BHs (17 ± 6 for — 150 GeV 
with a NPD spectrum) can be detected with GLAST. In this case, the DM signature will be the detection of up to 
20 point-like, bright, gamma-ray sources in a 3° circle around the Andromeda center (see Fig. 9). 



C. Contribution of SMBHs and IMBHs to the Gamma-ray Background 

EG RET data have shown the existence of an Extragalactive diffuse Gamma-ray Backgr ound (EGB) in the GeV 
range [lOOl Il0l| , initially interpreted as photons from an unresolved population of blazars |l02l | . Today such astro- 
physical sources are believed to play a significant role but they can hardly account for the entire background, while 
DM spikes an d mi ni-spikes could in principle strongly contribute to the fiux. 

UUio et al. [l03j estimated the contribution of DM annihilations from halos without central BHs to the EGB. The 
total number of photons that contributes to the EGB from redshift z is in this case: 

dM, - j dMH + z)3_(M, .)^^ • (29) 

Jdrp^ir,M,z)e-^i-o.)^^,E^dto, 

where dn/dM is the comov ing number density of DM halos with mass M at redshift z and it is computed using the 
Press-Schechter formalism [sol . Il04l |. The (1 + z)^ factor converts comoving to proper density, while dN.y/dE is the 
usual annihilation spectrum described in the last Section, p is the DM density and the factor e~'^'^^°'^^ takes into 
account the gamma ray absorption due to pair production with the ExtraGalactic Background Light in the infrared 
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Figure 9: Map of the gamma-ray flux, in units of photons cm~'^s~^, from DM annihilations around IMBHs in M31, relative to 
one random realization of IMBHs in M31. The size of the bins is 0.1° and the energy threshold is 4 GeV, as appropriate for 
GLAST. The big circle shows, for comparison, the M31 scale radius r^, and the small circles highlight IMBHs within the reach 
of GLAST (5(T detection in 2 months). Figure taken from Fornasa, Taoso and Bertone (2007). 



or optical band (r « (£^o/10 GeV)°-^(z/3.3)). is the energy today so that dEdt = d£'o(l + z)rfto(l + and Rq 
is the today scale factor. 

Each halo is described by a NFW (Eq. [1]) or M99 ( Eg. p ) profile, whose concentration parameter c {c — Tyir/i^s) 
depends on redshift as in Eke, Navarro and Steinmetz [lOSj . In order to reproduce the EGRET diffuse flux, a large 
boost factor would be required, that, if applied to the Galactic center, would grossly exceed the EGRET measurement 
of the gamma ray flux from Sgr A* [80] . 

The presence of a spike would only marginally affect the signal from the Galactic center, because even if the spike 
was present in the past for the SMBH of the MW, it would already be destroyed. However, for the EGB spikes are 
indeed important, since the background receives contributions from halos at hi gh re dshift, at a time when astrophysical 
and particle physics effects did not have time yet to damp DM overdensities [78]. 

Ahn et al. [78] included spikes around SMBHs in the computation made by UUio et ai. T he a uthors used the 
phenomelogical relation between the mass of the BH M, and the mass of the hosting halo Mh |l06l |: M./IO^M© = 
a{Mh/10^^ Mq)^ , where M/j is considered at z = and different values for (a, b) are used [78*1. 

For each halo at z = 0, a SMBH with mass Af, (related to M^) is posed at the redshift of BH formation z = z, and 
the halo is evolved till now (assuming that Mh remains unchanged). The spikes form and then decay in a self-similar 
way (see Section llV C[ p{r, t) « p{r, 0)k — p{r, 0)e~'^/^), so that the evolution of the spike is described by 

Psp{r,t) = Pspfin" (—^ , (30) 

with 7 and are the slopes of the initial and final profile, e ~ icisp — 7)"^ and r is the time since spike formation 
in units of the heating time theat- 

If we specify a NFW profile for the DM density before the spike formation, the complete DM distribution can be 
written as follows: 

pir)NFW {z > 2.) 

p{r)NFW (z < z,,r > rs,o) ■ (32) 

^ p{r)NFW + Psp{r, t) « Psp{r, t) {z < z„r < r^fi) 
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Figure 10: Left: Gamma-ray background produced by DM annihilations in DM halos with spikes (solid lines), compared to the 
halo-only contribution (dotted). The EGRET diffuse flux limits are shown for comparison. Three choices for the parameters 
in the relation between tha mass of the BH and the mass of the halo are used: {a,b) — (0.027, 1.82) in red, (0.10, 1.65) in 
blue and (0.67, 1.82) in green. The DM parameters adopted are m-^ = 100 GeV, av — 10^^''cm''s^^. Figure taken from Ahn, 
Bertone, Merritt and Zhang (2007) Right: Contribution to the EGB from scenario A and scenario B IMBH mini-spikes. Also 
shown are the EGRET data and predictions for the case of standardhalos without spikes and mini-spikes. The la scatter in 
the number of IMBHs is shown for scenario A. For scenario B, the la scatter for the number of IMBHs and for the upper limit 
in the redshift integration are shown combined. Results are shown for — 100 GeV and av — 3 ■ 10~^''cm^s~^. Figure taken 
from Horiuchi and Ando (2006). 



For a DM particle of 100 GeV the contribution to the EGB is summarized in the left panel of Fig. 10, that 
shows an enhancement of more than an order of magnitude with respect to the case without SMBHs. The main 
contribution comes from low-energies, due to the annihilations at high redshift in the still-present spikes. The energy 
of the produced photons is then redshifted till now. 

The EGB from DM annihilation in mini-spikes around IMBHs has been calculated in Ref. [t^. The results are 
shown in the right panel of Fig. 10 for a neutralino of 100 GeV and both the formation scenarios. As before, the 
main contribution is at low energies, but now the predictions from scenario B are able to account for all the EGB, 
and the same is true if one consider the line spectrum. This model is sensitive to the average of the halos properties, 
whereas, in the case of annihilation e.g. from the Galactic Center, one has to deal with a single realization that may 
differ significantly from the average, so that the study of gamma ray background and the way DM contributes is, for 
Indirect searches, an interesting alternative to the study of unidentified sources. 



VI. CONCLUSIONS 



We have reviewed the consequences of the growth and evolution of BHs on the distribution of stars and DM around 
them. Initially, the BH mass was considered to be fixed in time, with particles evolving around it under the influence 
of its gravitational potential and their gravitational collisions. A steady-state solution is predicted for stars only for 
collisional nuclei: it is a zero-flux solution characterized by a central cusp with a slope of —7/4 [i^ within 0.2 rh. 
DM particles, being coUisionless, are not affected by self gravitational interactions but by encounters with stars l47|: 
in the same time that the —7/4 cusp forms for stars, the DM distribution is rearranged in a profile with a central 
slope of —3/2, that will decay due to heating. Other effects can significantly reduce collisional overdensities, e.g. both 
stellar and DM cusps react to the presence of a BH binary [H^l and the stellar number density can be influenced by 
loss-cone dynamics 2^1. 

In the case of adiabat ic g rowth of a BH, instead, the final DM and stellar overdensities are much steeper, reaching 
slope as large as —2.25 j7ll. It^. But again, the evolution of the host galaxy drives a suppression of the spike, due to 
cumulative mergers [s^ and possible off-center BH formation [t^ . Even if DM overdensities around BHs may seem 
to be very promising for Indirect detection, the presence of spikes is probably not realistic (TSj . Alternative strategies 
include the study of DM contribution to the EGB and the search for IMBHs. In the first case, the gamma-ray 
background receives contribution from BHs at high redshift, when spikes and mini-spikes were already present, but 
damping mechanisms were not effective yet. DM annihilations from spikes around SMBHs at high redshift contribute 
significantly to the EGB [78|. A population of IMBHs is then able to reproduce the value detected by EGRET, even 
for standard assumptions for the initial DM profile (79| . 
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Figure 11: Exclusion plot in the DM mass vs. annihilation cross section plane. Above the thick (dotted) solid line all the 
statistical realizations of a MW-like halo in Bertone, Zentner and Silk (2005) contain at least one IMBH with an annihilation flux 
larger than the estimated point source sensitivity of GLAST (EGRET) in 2 months. Above the thin solid line all the statistical 
realizations of the Andromeda galaxy in Fornasa, Taoso and Bertone (2007) have at least one IMBH with an annihilation flux 
larger than 10~^''cm~'^s~^. The dotted line indicates the region above which the contribution from DM annihilation around 
IMBHs at high redshift can account for the EGB as measured by EGRET. The empty (full) dots represent SUSY models 
(obtained with DarkSusy) where the neutralino relic density is within 3 {5)a from the WMAP-measured value. 

Mini-spikes around IMBHs can effectively be searched for with GLAST and can be identified with a careful spec- 
tral analysis of the signal [13] ■ They can also been searched for in the Andromeda galaxy, distributed in a quite 
characteristic fashion that leaves no room for alternative astrophysical interpretations [9q . 

We show in Fig. 11 the constraining power of the different detection strategies discussed above: 

• the dotted and solid lines are relative to the population of IMBHs in our Galaxy [l^] ■ For each value of the mass 
of the DM candidate, we estimate the value of av for which all the aformentioned realizations of the Milky Way 
have at least one IMBH that should have been detected by EGRET at 5cr, as discussed in Ref. [l^l- In other 
words, if none of the EGRET unidentified sources can be interpreted as IMBH, a model with a combination 
{m^,av) above the dotted line is compatible with the data with a probability less than 0.5%. The solid line 
shows the prospectesd for detection with GLAST 

• the thin solid line is relative to the population of IMBHs in M31 [93]. For each value of we determined the 
value of av for which all the 200 realizations of M31 have at least one IMBH with an annihilation flux larger 
than the 5a sensitivity of GLAST (parametrized as in Ref. [95j). for an exposure of 2 months, and an energy 
threshold of 4 GeV. Null searches of IMBHs in Andromeda will exclude the region above the thin line 

• the dashed line is relative to the contribution of IMBHs to the EGB [t^ . Each point above the line corresponds 
to a DM candidate for which anni hilations in IMBHs mini-spikes at high redshift exceeds the EGB values as 
measured by EGRET in Ref. [lOOt . 

• the empty dots represent Supersymmetric Models (generated with DarkSusy jl07l |l for which the neutralino can 
account for all the DM, being within 3a from the WMAP determination of the DM relic density jl08l | 



• the full dots represent SUSY Models within 5a from WMAP 
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GLAST can thus probe regions of the parameter space at very low cross section, so that IA/[BHs can be soon 
discovered or ruled out. If none of these sources is detected by GLAST, it is very unlikely than such scenario will 
survive, at least for a WIMP DM candidate. 

The gam ma -ray flu x fro m DM annihilations in spikes and mini-spikes has been widely discussed in literature (Refs. 
[47|.150|. ItH. Wa. [tgI . [i09| . just to name a few) and many experimental collaborations are actively working to improve 
their discovery potential. But so far, the gravitational effects of BHs have only been observed through their impact on 
stellar populations, since there is not a detected signal currently and clearly interpretated as DM annihilation. Photons 
from DM annihilations could have been detected in the past as subdominant contributions from gamma-ray sources 
PH but the main challenge remains to disentagle such a contribution from other photon production mechanisms. 

Instead, the effect of a central BH on the stellar population has be en detect ed for the MW [H, [4^ and such 
interpretation is proposed also for other galaxies of the Local Group [llOl encouraging further attempts of 

detecting BH effects on DM and, eventually, of detecting DM indirectly through its annihilations products. 
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